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Abstract 

As RHIC is entering the precision measurement era and the LHC is producing a copious amount 
of new data, the role of 3+1D event-by-event viscous hydrodynamics is more important than 
ever to understand the bulk data as well as providing the background for hard probes. For more 
meaningful comparison with the experimental data, it is also important that hydrodynamics be 
coupled to the hadronic afterburner. In this proceeding we report on preliminary results of cou- 
pling MUSIC with UrQMD. 



1. Introduction 

Full understanding of Quark-Gluon Plasma (QGP) in relativistic heavy ion collisions requires 
a wide range of theoretical tool sets. There is little doubt that the initial condition of the nuclei 
before the collision as well as the gluon dynamics right after the collision is governed by clas- 
sical Yang-Mills dynamics JTJ [2j . Less than 1 fm/c after the collision, the system reaches the 
local equilibrium stage and hydrodynamic expansion begins. How the system reaches the local 
equilibrium state so quickly is still hotly debated, but progress is being made [3|. In the initial 
stage of the hydrodynamic expansion, the system is in the QGP phase. This phase is so dense that 
hydrodynamics must be valid. As the system expands and cools, it enters the hadronic phase. On 
the one hand, this system should still admit a hydrodynamic description via the hadronic equation 
of state. On the other hand, as the system becomes more dilute, hadronic kinetic theory models 
should become more accurate. Therefore, there should be a range of temperatures where both 
descriptions are valid. This enables us to switch from one description to another within a range 
of temperature slightly below the critical temperature. In this short proceeding, we report on 
our effort to deal with this transition between the hydrodynamics phase and the hadronic kinetic 
theory phase by coupling 3+ ID event-by-event viscous hydrodynamics (MUSIC) and a hadronic 
quantum molecular dynamics model (UrQMD). The initial Glasma phase and its coupling to the 
hydrodynamic phase is dealt with in Refs.(TJ|2). 

2. 3+1D Event-by-Event Viscous Hydrodynamics 

Hydrodynamics was first applied to particle collisions by Belenkij and Landau in 1956 [4|. 
In more recent times, boost-invariant 2+1D ideal hydrodynamics have enjoyed much success in 
heavy ion phenomenology OS- Detailed data from RHIC experiments, however, eventually 
made it clear that it is necessary to go beyond the 2+ ID ideal hydrodynamic with smooth av- 
eraged initial conditions. For instance, the surprisingly large triangular flow QE]|9l can never 
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Figure 1 : The momentum distribution of RHIC and the LHC 5 % most central collisions compared to the 
MUSIC+UrQMD calculations. RHIC result is based on 10,000 events (100 UrQMD events on each of 100 hydro events) 
events and LHC result is based on 1,000 events (10 UrQMD events on each of 100 hydro events). For RHIC, the results 
of pure MUSIC events are also shown. 




Figure 2: The pj dependent vi from RHIC and the LHC in the 10 - 20 % centrality class compared to the MUSIC and 
MUSIC+UrQMD calculations. RHIC result is based on 10,000 events (100 UrQMD events on each of 100 hydro events) 
events and LHC result is based on 1,000 events (10 UrQMD events on each of 100 hydro events). 



be reproduced in this way. More recent efforts have improved this by including non-zero shear 
viscosity OOlQj], 3+1D dynamics [El El, and fluctuating initial conditions OLE] ED. How " 
ever, none of the above mentioned studies implements all three improvements over the 2+ ID 
ideal hydrodynamics. Only recently, the first implementation of 3+1D event-by-event viscous 
hydrodynamics was made by 3 of the present authors ifTTl [181 [191 . 

Our implementation of hydrodynamics (named MUSIC) uses the hyperbolic r-r\ coordinate 
system without assuming invariance in r\. The viscous tensor if is calculated by using a variant 
of the Israel-Stewart formalism ||20l . The algorithm chosen to solve the system of hyperbolic 
equations is the Kurganov and Tadmor algorithm in the semi-discrete form ED together with 
Heun's method. This implementation is relatively simple yet very stable. Glauber approximation 
provides fluctuating initial conditions for the energy density. For the equation of state, the lattice 
based parameterization "s95p-vl" from [22 1 is the default option although other options are also 
available. When freezing out at about T = 140 MeV without the afterburner, the Cooper-Frye 
formula is used to get the final state hadronic spectrum with the full resonance decays adapted 
from Heinz and Kolb's 2-D scheme. This Cooper-Frye procedure includes the viscous correction 



triangular flow v3 (Air+Ari ZOOGeV, 1 0-20%) 
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Figure 3: The pj dependent triangular flow from RHIC and the LHC in the 10 - 20 % centrality class compared to the 
MUSIC and MUSIC +UrQMD calculations. RHIC result is based on 10,000 events (100 UrQMD events on each of 100 
hydro events) events and LHC result is based on 1,000 events (10 UrQMD events on each of 100 hydro events). The data 
points are from Refs.Q (PHENIX) and (9) (ALICE). 



3. UrQMD Afterburner 

UrQMD (Ultrarelativistic Quantum Molecular Dynamics) is a hadronic kinetic theory model 
that has been very successful in describing heavy ion collisions up to the RHIC energy [24, 
|251 . Recently, H. Petersen et al has updated the UrQMD code to include hydrodynamics [26 1 
and made it publicly available. The hydrodynamics included in the original version (UrQMD 
v. 3.3pl) is a 3+1D ideal hydrodynamic simulation in the Cartesian coordinate system with the 
simple freeze-out surface determined at a constant t. MUSIC on the other hand uses the hyper- 
bolic coordinate system and generates a fairly complicated iso-thermal freeze-out surface. We 
have implemented an interface between these two systems fully taking into account the freeze- 
out hypersurface directions. (See also Ref. ETll for an exact freeze-out-surface finding algorithm.) 
One missing ingredient is the viscous correction of the thermal distribution when it is sampled 
for the subsequent UrQMD run. This will be implemented in the near future although the viscous 
correction is not expected to be large at around T — 170 MeV for moderate values of pj ll23l . In 
this work, we present our preliminary results with previously tuned values of the hydrodynamics 
parameters. More complete analysis including hydro parameter dependencies will be reported 
elsewhere. 



4. Results and Summary 



Three main results are shown in Figures T|2 and[3] For RHIC, we have 10,000 events in each 



centrality bin generated by running 100 UrQMD events on each of 100 viscous MUSIC events. 
For the LHC, we have 1,000 events in each centrality bin generated by 10 UrQMD events on 
each of 100 viscous MUSIC events. The transition to UrQMD was made at T tv = 170 MeV. The 
two panels in Figure[T]show proof of principle results. The RHIC pj spectra for n~ ', K and p are 
reasonably well produced for the - 5 % centrality bin. For all 3 species, the description of the 
data is either improved over or of about the same quality as the pure hydrodynamics results. In 
the second panel, our calculations are compared with the ALICE data shown in this conference 
1281 . It should be emphasized here that the calculation was done before the data were shown. 

For V2(pr) in Figure [2] the p result is much improved with the afterburner at RHIC. This 
may be due to the better description of finite baryon mean free path in UrQMD compared to 
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hydrodynamics. For the LHC, one may argue that the afterburner improves the description a 
little, but the statistical error at this point is too big to say that definitely. Calculations with 
afterburner lead to larger v?,(pt) in general. However, this needs more careful study to quantify. 

In summary, we have reported our preliminary results on coupling MUSIC to the UrQMD 
afterburner. Encouragingly, changes from the pure hydrodynamics cases seem to go in the right 
direction. In near future, we plan to (i) improve our statistics up to 10 times, (ii) include the effect 
of 6f at the transition, and ( Hi) slightly re-tune hydrodynamics parameters such as the viscosity 
for better experimental fit. These further studies will be reported elsewhere. 

Acknowledgments 

CG, SJ and SR are supported by the Natural Sciences and Engineering Research Council of 
Canada. BPS is supported under DOE Contract No.DE-AC02-98CH 10886 and acknowledges 
support from a BNL Lab Directed Research and Development grant. BPS gratefully acknowl- 
edges a Goldhaber Distinguished Fellowship from Brookhaven Science Associates. We greatly 
appreciate computer time on the Guillimin cluster at the CLUMEQ HPC centre, a part of Com- 
pute Canada HPC facilities where bulk of these calculations were carried out. We also thank 
H. Petersen and P. Huovinen for helpful discussions. 

References 
References 

[1] C. Gale, S. Jeon, B. Schenke, P. Tribedy, R. Venugopalan, in: this proceedings. 
[2] C. Gale, S. Jeon, B. Schenke, P. Tribedy, R. Venugopalan (2012). arXiv : 1209 . 6330 
[3] A. Kurkela, G. D. Moore, JHEP 1 1 12 (201 1) 044. |arXiv: lToT. 5050 
[4] S. Belenkij, L. Landau, Nuovo Cim.Suppl. 3S10 (1956) 15. 

[5] P. F. Kolb , U. W. Heinz, in: R. C. Hwa, X.-N. Wang (Eds.), Quark Gluon Plasma 3, World Scientific, Singapore, 

2004. |arXiv : micl :: th /0305084 
[6] P. Huovinen, in: R. C. Hwa, X.-N. Wang (Eds.), Quark Gluon Plasma 3, World Scientific, Singapore, 2004. 

larXiv : mic l-th/0305064 
[7] E. Richardson (PHENIX Collaboration) (2012). arXiv : 1206 . 3501 
[8] Y. Pandit (STAR Collaboration) (2012). [arXlv: 1209 ■ 0244| 

[9] B. Abelev, et al. (ALICE Collaboration) (2012). |arXiv : 1205 ■ 5761] 

[10] U. W. Heinz, H. Song, A. K. Chaudhuri, Phys.Rev. C73 (2006) 034 904,|arXiv:nucl-th/0510014| 
[1 1] R. Baier, P. Romatschke, Eur.Phy s.J. C51 (2007) 677-687.|arXiv:micl-t h/0610108 
[12] S. Muroya, C. Nonaka (1997). [arXlv : nucl-th/9709004| 

[13] T. Hirano, Phys.Rev.Lett. 86 (2001) 2754-2757. |arXiv : nucl-th/0004029| 

[14] C. Aguiar, Y. Hama, T. Kodama, T. Osada, Nucl.Phys. A698 (2002) 63 9-642.|arXiv:hep-ph701 06266 
[15] H. Holopainen, H. Niemi, K. J. Eskola, Phys.Rev. C83 (2011) 0349 01. |arXiv: 1 007 . 0368 
[16] L. Pang, Q. Wang, X.-N. Wang, Phys.Rev. C86 (2012) 024911. arXiv: 1205.5019 
[17] B. Schenke, S. Jeon, C. Gale, Phys.Rev. C82 (2010) 014903. arXiv : 1004 . 1408 
[18] B. Schenke, S. Jeon, C. Gale, Phys.Rev.Lett. 106 (2011)042301. arXiv : 1009 . 3244 

[19] B. Schenke, S. Jeon, C. Gale, Phys.Lett. B702 (2011) 59-63. arXiv : 1102 . 0575 

[20] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, M. A. Stephanov, JHEP 0804 (2008) 100. |arXiv: 0712 .2451| 
[21] A. Kurganov, E. Tadmor, Journal of Computational Physics 160 (2000) 241. 

[22] P. Huovinen, P. Petreczky, Nucl.Phys. A837 (2010) 26-53. [arXiv : 0912 ■ 254l] 

[23] M. Dion, J.-F. Paquet. B. Schenke, C. Young, S. Jeon, et al., Phys.Rev. C84 (201 1) 064901 . [arXiv : 1 109 ■ 44"05| 
[24] S. Bass, M. Belkacem, M. Bleicher, M. Brandstetter, L. Bravina, et al., Prog.Part.Nucl.Phys. 41 (1998) 255-369. 
larXiv : mic l-th/9803035 

[25] M. Bleicher, E. Z abrodin, C. Spieles, S. Bass, C. Ernst, et al., J.Phys. G25 (1999) 1859-1896. 
arXi v : hep-ph/9909407 

[26] H. Petersen, J. Steinheimer, G. Burau, M. Bleicher, H. Stocker, Phys.Rev. C78 (2008) 044901 . [arXiv : 0806 ■ 1695] 

4 



[27] P. Huovinen, H. Petersen (2012). arXiv : 1206 . 3371 
[28] L. Milano (ALICE Collaboration), in: this proceedings. 



5 



